Study Profile of the Tsuruoka Metabolomics Cohort Study (TMCS).

The Tsuruoka Metabolomics Cohort Study (TMCS) is an ongoing population-based cohort study being conducted in the rural area of Yamagata Prefecture, Japan. This study aimed to enhance the precision prevention of multi-factorial, complex diseases, including non-communicable and aging-associated diseases, by improving risk stratification and prediction measures. At baseline, 11,002 participants aged 35-74 years were recruited in Tsuruoka City, Yamagata Prefecture, Japan, between 2012 and 2015, with an ongoing follow-up survey. Participants underwent various measurements, examinations, tests, and questionnaires on their health, lifestyle, and social factors. This study used an integrative approach with deep molecular profiling to identify potential biomarkers linked to phenotypes that underpin disease pathophysiology and provide better mechanistic insights into social health determinants. The TMCS incorporates multi-omics data, including genetic and metabolomic analyses of 10,933 participants and comprehensive data collection ranging from physical, psychological, behavioral, and social to biological data. The metabolome is used as a phenotypic probe because it is sensitive to changes in physiological and external conditions. The TMCS focuses on collecting outcomes for cardiovascular disease, cancer incidence and mortality, disability, functional decline due to aging and disease sequelae, and the variation in health status within the body represented by omics analysis that lies between exposure and disease. It contains several sub-studies on aging, heated tobacco products, and women's health. This study is notable for its robust design, high participation rate (89%), and long-term repeated surveys. Moreover, it contributes to precision prevention in Japan and East Asia as a well-established multi-omics platform.

East Asian-specific and cross-ancestry genome-wide meta-analyses provide mechanistic insights into peptic ulcer disease.

Peptic ulcer disease (PUD) refers to acid-induced injury of the digestive tract, occurring mainly in the stomach (gastric ulcer (GU)) or duodenum (duodenal ulcer (DU)). In the present study, we conducted a large-scale, cross-ancestry meta-analysis of PUD combining genome-wide association studies with Japanese and European studies (52,032 cases and 905,344 controls), and discovered 25 new loci highly concordant across ancestries. An examination of GU and DU genetic architecture demonstrated that GUs shared the same risk loci as DUs, although with smaller genetic effect sizes and higher polygenicity than DUs, indicating higher heterogeneity of GUs. Helicobacter pylori (HP)-stratified analysis found an HP-related host genetic locus. Integrative analyses using bulk and single-cell transcriptome profiles highlighted the genetic factors of PUD being enriched in the highly expressed genes in stomach tissues, especially in somatostatin-producing D cells. Our results provide genetic evidence that gastrointestinal cell differentiations and hormone regulations are critical in PUD etiology.

Integrated analysis of human DNA methylation, gene expression, and genomic variation in iMETHYL database using kernel tensor decomposition-based unsupervised feature extraction.

Integrating gene expression, DNA methylation, and genomic variants simultaneously without location coincidence (i.e., irrespective of distance from each other) or pairwise coincidence (i.e., direct identification of triplets of gene expression, DNA methylation, and genomic variants, and not integration of pairwise coincidences) is difficult. In this study, we integrated gene expression, DNA methylation, and genome variants from the iMETHYL database using the recently proposed kernel tensor decomposition-based unsupervised feature extraction method with limited computational resources (i.e., short CPU time and small memory requirements). Our methods do not require prior knowledge of the subjects because they are fully unsupervised in that unsupervised tensor decomposition is used. The selected genes and genomic variants were significantly targeted by transcription factors that were biologically enriched in KEGG pathway terms as well as in the intra-related regulatory network. The proposed method is promising for integrated analyses of gene expression, methylation, and genomic variants with limited computational resources.

Comparison of the loci associated with HbA1c and blood glucose levels identified by a genome-wide association study in the Japanese population.

Aims: Hemoglobin A1c (HbA1c) levels are widely employed to diagnose diabetes. However, estimates of the heritability of HbA1c and glucose levels are different. Therefore, we explored HbA1c- and blood glucose-associated loci in a non-diabetic Japanese population. Methods: We conducted a two-stage genome-wide association study (GWAS) on variants associated with HbA1c and blood glucose levels in a Japanese population. In the initial stage, data of 4911 participants of the Japan Multi-Institutional Collaborative Cohort (J-MICC) were subjected to discovery analysis. In the second stage, two datasets from the Tohoku Medical Megabank project, with 8175 and 40,519 participants, were used for the replication study. Association of the imputed variants with HbA1c and blood glucose levels was determined via linear regression analyses adjusted for age, sex, body mass index (BMI), smoking, and genetic principal components (PC1-PC10). Moreover, we performed a BMI-stratified GWAS on HbA1c levels in the J-MICC. The discovery analysis and BMI-stratified GWAS results were validated with re-analyses of normalized HbA1c levels adjusted for site in addition to the above, and blood glucose adjusted for fasting time as an additional covariate. Results: Genetic variants associated with HbA1c levels were identified in KCNQ1 and TMC6. None of the genetic variants associated with blood glucose levels in the discovery analysis were replicated. Association of rs2299620 in KCNQ1 with HbA1c levels showed heterogeneity between individuals with BMI ≥ 25 kg/m2 and BMI < 25 kg/m2. Conclusions: The variant rs2299620 in KCNQ1 might affect HbA1c levels differentially based on BMI grouping in the Japanese population. Supplementary Information: The online version contains supplementary material available at 10.1007/s13340-023-00618-0.

Epigenetic profile of Japanese supercentenarians: a cross-sectional study.

BACKGROUND: Centenarians and supercentenarians with exceptional longevity are excellent models for research towards improvements of healthy life expectancy. Extensive research regarding the maintenance and reduction of epigenetic age has provided insights into increasing healthy longevity. To this end, we explored the epigenetic signatures reflecting hallmarks of exceptional healthy longevity, including avoidance of age-related diseases and cognitive functional decline. METHODS: In this cross-sectional study, we enrolled Japanese non-centenarians (eligible participants aged 20-80 years) from the Tohoku Medical Megabank Community-Based Cohort Study and centenarians and supercentenarians (aged 101-115 years) from the Tokyo Centenarian Study and the Japanese Semi-supercentenarian Study. We assessed participants' whole-blood DNA methylation profiles and then developed sex-specific and non-specific first-generation epigenetic clocks by elastic net regression, calculated individuals' epigenetic ages, and assessed their age acceleration. We also screened for age-related CpG sites in non-centenarians by epigenome-wide linear regression analyses and ANOVA. We subsequently investigated which CpG sites in centenarians and supercentenarians had DNA methylation patterns following the age-related findings obtained from non-centenarians and which did not. We further characterised CpG sites with hypermethylation or hypomethylation in the centenarians and supercentenarians using enrichment and protein-protein interaction network analyses. FINDINGS: We enrolled 421 non-centenarians (231 [55%] women and 190 [45%] men; age range 20-78 years), recruited between May 20, 2013, and March 31, 2016, and 94 centenarians and supercentenarians (66 women [70%] and 28 [30%] men; age range 101-115 years), recruited between Jan 20, 2001, and April 17, 2018. Non-sex-specific epigenetic clock showed the highest accuracy (r=0·96) based on which centenarians and supercentenarians had negative epigenetic age acceleration. Epigenome-wide association analyses further showed that centenarians and supercentenarians had younger-than-expected epigenetic states (DNA methylation profiles similar to those of non-centenarians) for 557 CpG sites enriched in cancer-related and neuropsychiatric-related genes, whereas these individuals had advanced (or older) epigenetic states for 163 CpG sites represented by genes related to TGF-ß signalling, which is involved in anti-inflammatory responses and known to contribute to healthy ageing. INTERPRETATION: These results indicate that exceptionally healthy longevity depends not only on maintaining young epigenetic states but also on advanced states of specific epigenetic regions. FUNDING: The Japan Agency for Medical Research and Development, KDDI Research, and Keio University. TRANSLATION: For the Japanese translation of the abstract see Supplementary Materials section.

Development and validation of genome-wide polygenic risk scores for predicting breast cancer incidence in Japanese females: a population-based case-cohort study.

PURPOSE: This study aimed to develop an ancestry-specific polygenic risk scores (PRSs) for the prediction of breast cancer events in Japanese females and validate it in a longitudinal cohort study. METHODS: Using publicly available summary statistics of female breast cancer genome-wide association study (GWAS) of Japanese and European ancestries, we, respectively, developed 31 candidate genome-wide PRSs using pruning and thresholding (P + T) and LDpred methods with varying parameters. Among the candidate PRS models, the best model was selected using a case-cohort dataset (63 breast cancer cases and 2213 sub-cohorts of Japanese females during a median follow-up of 11.9 years) according to the maximal predictive ability by Harrell's C-statistics. The best-performing PRS for each derivation GWAS was evaluated in another independent case-cohort dataset (260 breast cancer cases and 7845 sub-cohorts of Japanese females during a median follow-up of 16.9 years). RESULTS: For the best PRS model involving 46,861 single nucleotide polymorphisms (SNPs; P + T method with PT = 0.05 and R2 = 0.2) derived from Japanese-ancestry GWAS, the Harrell's C-statistic was 0.598 ± 0.018 in the evaluation dataset. The age-adjusted hazard ratio for breast cancer in females with the highest PRS quintile compared with those in the lowest PRS quintile was 2.47 (95% confidence intervals, 1.64-3.70). The PRS constructed using Japanese-ancestry GWAS demonstrated better predictive performance for breast cancer in Japanese females than that using European-ancestry GWAS (Harrell's C-statistics 0.598 versus 0.586). CONCLUSION: This study developed a breast cancer PRS for Japanese females and demonstrated the usefulness of the PRS for breast cancer risk stratification.

Evaluation of short-term epigenetic age fluctuation.

Considerable effort has been spent on lowering and maintaining the epigenetic age. However, the extent to which epigenetic age fluctuates under normal conditions is poorly understood. Therefore, we analyzed methylation data from monocytes and peripheral blood mononuclear cells collected from two Japanese men. The ranges of the Pan-tissue, Skin and blood, and DNAm PhenoAge epigenetic age during 3 months were ≥ 5.62, ≥ 3.04, and ≥ 8.23 years, and the maximum daily changes were 5.21, 3.20, and 6.53 years, respectively. These fluctuations were not suppressed by correcting for cell-type composition. Although the underlying biological mechanism remains unclear, there was a nonnegligible degree of age fluctuation which should inform personalized clinical applications.

Low MICA Gene Expression Confers an Increased Risk of Graves' Disease: A Mendelian Randomization Study.

Background: Expression of natural killer group 2 member D (NKG2D) ligand (NKG2DL) plays a major role as a "danger signal" on stressed cells to promote removal of the latter by NKG2D-expressing cytotoxic lymphocytes. NKG2DL expression has been found in peripheral immune cells as well, such as in macrophages; however, the effect of this expression is yet to be determined. Methods: We determined instrumental variables (IVs; R2 <0.01 in linkage disequilibrium), explaining the major variance in major histocompatibility complex class I chain-related protein A (MICA) and B (MICB) gene expression levels from the expression-quantitative trait locus (eQTL) of NKG2DLs based on the RNA-seq analysis of peripheral blood mononuclear cells (PBMCs) from 381 Japanese. Simultaneously, the target outcomes were filtered by PheWAS from 58 health risks, using a community-based cohort study composed of 44,739 Japanese residents. Finally, we estimated the causal effect of gene expression levels on the outcomes using the Mendelian randomization approach. Results: We determined nine and four IVs, explaining 87.6% and 33.0% of MICA and MICB gene expression levels, respectively. In the association test, we identified 10 or 13 significant outcomes associated with the MICA or MICB eQTLs, respectively, as well as the causal effect of MICA expression on Graves' disease (GD) (p = 4.2 × 10-3; odds ratio per 1 S.D. difference in the expression: 0.983 [confidence interval: 0.971-0.995]), using the weighted median estimator, without significant pleiotropy (p > 0.05), and the results were consistent across the sensitivity analyses. Conclusions: Our study provide novel evidence associating NKG2DL expression with GD, an autoimmune thyroiditis; direction of the effect indicated the immunoregulatory role of MICA expression in PBMCs, suggesting the importance of further functional assays in inflammatory diseases.

Longitudinal DNA methylation dynamics as a practical indicator in clinical epigenetics.

BACKGROUND: One of the fundamental assumptions of DNA methylation in clinical epigenetics is that DNA methylation status can change over time with or without interplay with environmental and clinical conditions. However, little is known about how DNA methylation status changes over time under ordinary environmental and clinical conditions. In this study, we revisited the high frequency longitudinal DNA methylation data of two Japanese males (24 time-points within three months) and characterized the longitudinal dynamics. RESULTS: The results showed that the majority of CpGs on Illumina HumanMethylation450 BeadChip probe set were longitudinally stable over the time period of three months. Focusing on dynamic and stable CpGs extracted from datasets, dynamic CpGs were more likely to be reported as epigenome-wide association study (EWAS) markers of various traits, especially those of immune- and inflammatory-related traits; meanwhile, the stable CpGs were enriched in metabolism-related genes and were less likely to be EWAS markers, indicating that the stable CpGs are stable both in the short-term within individuals and under various environmental and clinical conditions. CONCLUSIONS: This study indicates that CpGs with different stabilities are involved in different functions and traits, and thus, they are potential indicators that can be applied for clinical epigenetic studies to outline underlying mechanisms.

Characterization of a novel nicotinamide adenine dinucleotide-cytochrome b5 reductase mutation associated with canine hereditary methemoglobinemia.

Hereditary methemoglobinemia associated with nicotinamide adenine dinucleotide-cytochrome b5 reductase (b5R) deficiency is a rare autosomal recessive disorder in animals. Recently, nonsynonymous b5R gene (CYB5R3) variants have been reported to be associated with canine and feline hereditary methemoglobinemia. However, the underlying molecular mechanisms of canine and feline methemoglobinemia caused by these nonsynonymous variants have not yet been reported. Previously, we reported a Pomeranian dog family with hereditary methemoglobinemia, carrying CYB5R3 mutation of an A>C transition at codon 194 in exon 7, replacing an isoleucine residue with leucine (p.Ile194Leu). In this study, we investigated the enzymatic and structural properties of the soluble form of wild-type and Ile194Leu canine b5Rs to characterize the effects of this missense mutation. Our results showed that the kinetic properties of the mutant enzyme were not affected by this amino acid substitution. The secondary structure of the wild-type and Ile194Leu b5Rs detected by circular dichroism showed a similar pattern. However, the mutant enzyme exhibited decreased heat stability and increased susceptibility to trypsin hydrolysis. Moreover, the thermostability and unfolding measurements indicated that the mutant enzyme was more sensitive to temperature-dependent denaturation than the wild-type b5R. We concluded from these results that unstable mutant enzyme properties with normal enzymatic activity would be associated with hereditary methemoglobinemia in the Pomeranian dog family.

Retrospective study on intercurrent pancreatitis with Babesia gibsoni infection in dogs.

Babesia rossi infection has been reported to be associated with the high prevalence of pancreatitis in dogs. In this study, we retrospectively investigated whether pancreatitis occurs in B. gibsoni-infected dogs. The clinical manifestations, and hematological and serum biochemical examination results, including canine pancreatic-specific lipase (cPL), in 20 B. gibsoni-infected dogs were analyzed. The cPL concentration exceeded 400 µg/l in only 2 dogs, and they were suspected of having pancreatitis. Although the cPL concentration did not correlate with the degree of anemia or the level of parasitemia, it correlated with the band neutrophil count, platelet count, and blood urea nitrogen (BUN) level. Our study suggested that the prevalence of pancreatitis is lower among B. gibsoni-infected dogs than B. rossi-infected dogs.